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ABSTRACT 

We evaluate f n /m p , fx/^rip, 1/fp, an d m^j (f^nip), extrapolated to 
physical quark mass, zero lattice spacing and infinite volume, for lattice 
QCD with Wilson quarks in the valence (quenched) approximation. The 
predicted ratios differ from experiment by amounts ranging from 12% to 
17% equivalent to between 0.9 and 2.8 times the corresponding statistical 
uncertainties. 
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1 INTRODUCTION 



In a recent paper [|IJ we presented lattice QCD predictions for the masses of 
eight low-lying hardrons, extrapolated to physical quark mass, zero lattice spacing, 
and infinite volume, using Wilson quarks in the valence (quenched) approximation. 
The masses we found were within 6% of experiment, and all differences between pre- 
diction and experiment were consistent with the calculation's statistical uncertainty. 
We argued that this result could be interpreted as quantitative confirmation of the 
low-lying mass predictions both of QCD and of the valence approximation. It ap- 
peared to us unlikely that eight different valence approximation masses would agree 
with experiment yet differ significantly from QCD's predictions including the full 
effect of quark- ant iquark vacuum polarization. 

We have now evaluated the infinite volume, zero lattice spacing, physical quark 
mass limit of f n /m p , fK/ m P , 1/fp, an d Tri^/{f^m^). To our knowledge there have 
been no previous systematic calculations of this physical limit of lattice meson decay 
constants. A review of earlier work is given in Ref. 0. The predicted m^j '(f^rrip) 
lies above its observed value by about its statistical uncertainty of approximately 
15%. The predicted f w /m p and 1/ f p are close to 12% below experiment, equivalent 
to about 0.9 times the corresponding statistical uncertainties, and fx/^p is below 
experiment by 17%, equal to 2.8 times its statistical uncertainty. 

Although overall the predicted values could be considered in fair agreement 
with experiment, the result that three of four decay constants range from 0.9 to 2.8 
standard deviations below experiment suggests the possiblity that the true, underlying 
values, if determined without statistical uncertainties, would fall somewhat below 
experiment. One possible source of such disagreement is that we have assumed exact 
isospin symmetry and ignored electromagnetic effects. In particular, although our 
prediction for f n lies below measured values determined from charged pion decays, 
reported numbers for the neutral pion decay constant |§ are quite close to our result. 
The significance of the closer agreement between our predicted f n and the neutral 
pion decay constant is unclear to us, however, as a consequence of the systematic 
uncertainties arising in the experimental determination of the neutral pion decay 
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constant. 

Another possible source of disagreement between our numbers and experiment 
may be the valence approximation itself. A simple physical argument tends to support 
of this alternative The valence approximation may be viewed as replacing the 
momentum and frequency dependent color dielectric constant arising from quark- 
antiquark vacuum polarization with its low momentum limit || . At low momentum, 
then, the effective quark charge appearing in the valence approximation will agree 
with the low momentum effective charge of the full theory. The valence approximation 
might thus be expected to be fairly reliable for low-lying baryon and meson masses, 
which are determined largely by the long distance behavior of the chromoelectric 
field. The valence approximation's effective quark charge at higher momentum can 
be obtained from the low momentum charge by the Callan-Symanzik equation. As a 
consequence of the absence of dynamical quark-antiquark vacuum polarization, the 
quark charge in the valence approximation will fall faster with momentum than it 
does in the full theory, at short distance the attractive quark-antiquark potential in 
the valence approximation will be weaker than in the full theory, and meson wave 
functions in the valence approximation will be pulled into the origin less than in the 
full theory. Since decay constants are proportional to the square of wave functions 
at the origin, decay constants in the valence approximation could be expected to be 
smaller than in the full theory. 

The calculations described here were done on the GF11 parallel computer at 
IBM Research and use the same collection of gauge configurations and quark 
propagators generated for the mass calculations of Ref . || . The full set of mass and 
decay constant calculations took approximately one year to complete. GF11 was used 
in configurations ranging from 384 to 480 processors, with sustained speeds ranging 
from 5 Gflops to 7 Gflops. With the present set of improved algorithms and 480 
processors, these calculations could be repeated in less than four months. 
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2 DEFINITIONS 

The normalization we adopt for pseudoscalar and vector decay constants in 
continuum QCD is 

<0\J?{0)\V(p,e,j)> = e»m 3 F 3) 

<0|jf(0)|P(j9,j)> = p"/* 

for vector and pseudoscalar states, \V(p, e, j) > and \P(p, j) >, respectively, normal- 
ized by 

<p\q> = (2rr) 3 p°5(p- q). 

Here j is a flavor-SU(3) octet index and vector and axial vector flavor-SU(3) currents 
Jj(x) and Jj^ix) are related to cannonical continuum quark and antiquark fields, ip c 
and ipc, and an orthonormal set of flavor-SU(3) matrices Xj by 

J?(x) = ^ c {x)^X^ c {x), 
Jf{x) = ^ c (x)7 5 7^A^ c (x), 

triXjXj) = -. 
Assuming exact isospin symmetry, we have 

fi = fn £ = 1, . . . 3, 

fi = fx i = 4, • • • 7, 
F l = ^ i = l,...3. 

In the valence approximation we have, in addition, 

8 V2U 

For simplicity, we will also use the names F p and for i<\ and F&, respectively. Our 
normalization gives f n the value (93.15 ± 0.11) MeV. 

The hadron mass calculation of Ref. Q was done using gaussian smeared 
lattice quark and antiquark fields defined in Coulomb gauge. Smeared fields have, 



3 



therefore, also been adopted for the present calculation. The smeared field <p r (x,t) is 
related by the lattice antiquark field if)g(x,t) by 

y 

I 7*1 2 

G r (z) = (\/irr)~ 3 exp(- 



The field <f> r (x, t) is defined correspondingly from ip £ (x, t). We take the smeared fields 
(f)Q (x) and (f) (x) to be ip(x) and if)(x), respectively. From these fields, define smeared 
currents by 

Jf r {x) = <j> r (x)'y s \j<f> r (x), 
Jf r {x) = r (x)7 /i A i r (x), 
J^(x) = ^(xhVA.Mx). 

Define the correlation functions C^ r , r (t), Cj r , r (t) and Cf r , r (t) to be 

cf r , r {t) = x:<[4,(f,t)]t4( ,o)>, (2.i) 

X 

Cj r , r {t) = X:<[Jj r ,(f,t)]tjj r (0,0)>, (2.2) 

X 

Cf, r (t) = ^<[J$(x,t)?J™(p,0)>. (2.3) 

X 

Then for a set of constants Z?, r , Zj r , T and Zf T , ry these correlation functions have the 
asymptotic forms, for a large time separation t and lattice time direction period T, 

Cf r , r {t) -> Z^., r {exp(—Tnft) + exp[—mf(T — t)]}, (2.4) 
C£ r (*) -> Z] / r , r {exp(-mJt) + exp[-m] / (T-t)]} (2.5) 
C£, r (t) -> ^; v {exp(-mft) + ea^[-mf(T-t)]}. (2.6) 

Here mj 3 and are pseudoscalar and vector masses, respectively. Eqs. fl2.4|) and 
(|2.6|) imply, in addition, the asymptotic form 



w j r / r //\ y L J Zj j r / r tt 



(2.7) 



jr'r 



Measured in units of the lattice spacing a the decay constants fja and Fja are 
then given, for any choice of smearing size r, by 

ni AyAP\2 

(fja? = (2.8) 

/ / by Hi Z—i j 

The coefficients and zj are finite renormalizations chosen so that the lattice cur- 
rents zfa 3 JjQ and zJa 3 J^ approach the continuum currents and Jj 1 , respectively, 
as the lattice spacing approaches zero. These constants are often given the "naive" 
values 

AN _ 9k P 
Zj - ZKj , 

zJ N = 2kJ. (2.10) 

where and kj are the hopping constants corresponding to the mass of the quark 
and antiquark for a pseudoscalar or vector meson, respectively, with flavor j, assumed 
here to have m q = m-q. A standard derivation of the naive finite renormalization 
constants, however, which relates the quark terms in the lattice action to those in the 
continuum action is actually not correct ||. Altough the naive finite renormalization 
constants do lead to the correct continuum limit, naive renormalization contributes 
to decay constant an additional extraneous dependence on lattice spacing. Numerical 
evidence for this behavior will be given in Section [F[ A consistent calculation of finite 
renormalizations ||, correct to zeroth in a mean- field theory improved perturbation 
expansion, gives 

Zj - I 



Ak r 



v 



4k, 



where k c is the critical hopping constant at which the pion's mass becomes zero. To 
first order in improved perturbation theory, the finite renormalizations become ZjZ 
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and zJ 1 z V0 where 

Z A1 = l-0.31cw(l/a), 

Z V1 = 1-0.820^(1/0). (2.12) 

Decay constants for mesons with m q ^ m-q will be discussed below. 

For equal values of the up and down quark masses, m u and m d , we adopt the 
convention that a missing flavor index j always has the value 1. Thus 

C£ P (t) = C[ r , r (t), (2.13) 

yP yP 

C r , r (t) = C lr , r (t), 

yV _ yV 

tyA _ tyA 



3 PROPAGATORS 

Table [I] lists the lattice sizes, parameter values, sweeps skipped between gauge 
configurations, and number of configurations used in the ensembles from which decay 
constants were calculated. Gauge configurations were generated with the Cabbibo- 
Marinari-Okawa algorithm. A variety of different test support the expectation that 
the large number of sweeps run between successive configurations were more than suf- 
ficient to produce statistically independent values for all of the quantities required by 
the present calculations. A discussion of the algorithms by which quark propagators 
were found is given in Ref. Q. 

Vector and pseudoscalar masses and the coefficients Zfi r , Z^, r , and Zft r were 
then determined by fitting Cfi r (t), Cy r (t), and Cft r (t) to the asymptotic forms of Eqs. 
(|2.4j) - ( p.6|) . As an alternative means of determining Zfi r given a value of Zfi r , fits 
of Cft r (t)/Cr' r (t) were also made to the asymptotic form of Eqs. Q2.7|). 
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For the lattice 8 3 x 32 at (3 of 5.70 we calculated propagators only for source r 
and sink r' of size 0. In all other cases we calculated propagators only for source size 
r of 2. To determine decay constants, according to Eqs. (p.8[) and (|2.9|) , fits for the 
lattice 8 3 x 32 to propagators for the single sink of size are sufficient, while for all 
other lattices fits are needed for sink sizes of both and 2. To determine the range of 
time separations to be used in fitting for each (3 and k value, we evaluated effective 
masses m p (t), m v (t) and m p ' (t) by fitting C P r (t), C^ r (t), and Cy r (t), respectively, 
to Eqs. ( |2.4|) - ( |2.6| ) at time separations t and t + 1. The largest interval at large 
t showing an approximate plateau in an effective mass graph we chose as the initial 
trial range on which to fit each propatator to the corresponding asymptotic form of 
Eqs. Q2.4|) - (|2.6| ). Similarly, the largest interval at large t showing an approximate 
plateau in a graph of Cft r „ (t) / C P r (t) we chose as the initial trial range for a fit to Eq. 
( |2.7| ). Figures ([[]) - (|1|) show the plateaus at large t in the effective masses m p (t), 
m p '(t), the plateau at large t in the ratio C^(t)/C p 2 (t), and the effective mass m v (t) 
for propagators with source r of 2 and sink r' of 0. Figures (D - (D show results 
for the lattice 16 3 x 32 at (3 of 5.70 and the largest corresponding k value, 0.1675. 
Figures (§) - (|8|) show results for the lattice 24 3 x 36 at (3 of 5.93 and the largest 
corresponding k, 0.1581. Figures (||) - (|12|) show results for 30 x 32 2 x 40 at (3 of 6.17 
and the largest corresponding k, 0.1532. 

Fits to data for a range of t were done by minimizing the fit's x 2 determined 
from the full correlation matrix for the data being fit. An automatic fitting program 
repeatedly carried out fits on every connected interval of four or more points within 
the initial trial fitting range. For fits to Eq. ( |2.7|) , which require only a single fitting 
parameter, we looked at intervals of three or more points. The final fitting range was 
chosen by the program to be the interval with the smallest value of x 2 P er degree 
of freedom. Altough a variety of other criteria could be used to determine the final 
fitting range, an advantage of the method we adopted is that it can be implemented 
automatically thereby reducing the potential for biases. The reliability of our in our 
final extrapolated results depends to some degree on adopting a consistent choice of 
fitting ranges at different parameter values. 

The horizontal lines in Figures (|1|) - ( |i~2"|) show the fitted values of masses and 
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Zq 2 /Zq 2 , and the pair of vertical lines in each figure indicates the interval of t values in 
the final fitting range chosen. It is perhaps useful to mention that since the effective 
mass at each t depends on data both at t and t + 1, the effective mass shown at the 
highest t within each final fitting range depends on data outside the fitting range. 
Thus the fitted lines tend to approximate the average of the effective masses within 
the fitting range but with the effective mass at highest t omitted. In all but one 
case, the data shows clear plateaus at large t extending over more than four time 
values. These plateaus appear to be fitted fairly reliably. For the rho propagator, 
C^(t), on the lattice 24 3 x 36, the plateau is more ambiguous. The fit is made over 
the four t values from 9 to 12, for which the three corresponding effective masses fall 
slowly. The effective mases at t of 12 to 15 then fall about one standard deviation 
below the fit. At t of 16 and 17, the effective masses then return to the original fitted 
value. This behavior is consistent with a statistical fluctuation although it makes the 
identication of the plateau at which to fit more ambiguous. Comparable ambiguities 
do not occur elsewhere in our data. The interval chosen is approximately a rescaling 
of the rho fitting intervals of 5 to 8, at (3 of 5.70 and 13 to 16, at (3 of 6.17. Also, the 
rho decay constant obtained from this fit in Section |] is interpolates smoothly values 
obtained from less ambiguous fits at other values of quark mass and lattice spacing. 
If this point were simply elimated from our extrapolations of f p in quark mass and 
in lattice spacing, our final continuum predictions would be nearly unchanged. 

Value of x 2 f° r t ne fits in Figures ([!]) - (|l^) are shown in Table [2|. Our fits 
for sinks with size r' of 2, fits at smaller k, and fits on the lattice 24 3 x 32 are of 
comparable quality to those shown and give comparable \ 2 . 

Statistical uncertainties of parameters obtained from fits and of any function of 
these parameters were determined by the bootstrap method O. From each ensemble 
of N gauge configurations, 100 bootstrap ensembles were generated. Each bootstrap 
ensemble consists of a set of N gauge configurations randomly selected from the 
underlying N member ensemble allowing repeats. For each bootstrap ensemble the 
entire fit was repeated, including a possibly new choice of the final fitting interval. 
The collection of 100 bootstrap ensembles thus yields a collection of 100 values of any 
fitted parameter or any function of any fitted parameter. The statistical uncertainty 
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of any parameter is taken to be half the difference between a value which is higher 
than all but 15.9% of the bootstrap values and a value which is lower than all but 
15.9% of the bootstrap values. In the limit of large N the collection of bootstrap 
values of a parameter p approaches a gaussian distribution and the definition we use 
for statistical uncertainty approaches the dispersion, d, given by y/< p 2 > — < p > 2 . 

In the absence of some independent method for determing the predictions of 
QCD, it appears inevitable that the choice of t interval on which to fit data to a large 
t asymptotic form must be made by some procedure which depends on the Monte 
Carlo data itself. Thus the statistical uncertainties in the data lead to a corresponding 
uncertainty in the choice of fitting interval which, in turn, leads to some additional 
uncertainty in the fitted result. Another advantage of our procedure for choosing the 
fitting interval combined with bootstrap evaluation of statistical uncertainties is that 
the values we obtain for statistical uncertainties include the uncertainty arising from 
the choice of fitting interval. A comparison of the error bars found for our final fits 
with the error bars found using the same fitting range held fixed across the bootstrap 
ensemble shows that typically about 10% of the final statistical uncertainty comes 
from fluctuations over the bootstrap ensemble of the fitting range itself. 



4 DECAY CONSTANTS 

To construct decay constants using the finite renomalizations of Eqs. ( |2.11| ) or 



fl2.12|) , we require values of the critical hopping constant k c and of the strong coupling 
constant a-msil/a). For k c we used the values determined in Ref. 0. These are listed 
in Table To determine a^l/a) we first calculated a-msi^ /a) using the mean-field 
improved perturbation theory relation || £| 

1 < TrU/3 > , N 

2 ( , s = t- + 0-025, 4.1 

9ws( 7r / a ) 9lat 

where < TrU > is the expectation value of the trace of a plaquette and gf at is 6/(3. We 
then found c%p{7r/a) given by Atx / gj^iji / a) and used the two- loop Callan-Symanzik 
equation to determine «ms(l/a). The corresponding values of z Al and z vi are shown 
in Table (HI). 
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Values of the decay constant in lattice units f w a for the various lattices, (3 and 
k shown in Table ([1]) are listed in Tables © - © ■ As a measure of the lattice spacing 
in each case, Table (0) gives the rho mass in lattice units m p (m n )a, extrapolated to 
the "normal" quark mass m n which produces the physical value of m n (m n )/m p (m n ) 
[p]]. We list also values of m n itself. These parameters are not given for the lattice 
8 3 x 32 since in this case we were not able to calculate propagators at small enough 
quark mass to perform the required extrapolation reliably. The finite renormalizations 
for the decay constants in Tables (|]) - (0) all include both the leading term z A0 of 
Eqs. (|2.11|) and the first order mean-field improved perturbative correction z A1 of 
Eqs. ( |2.12D . The second column in each table gives values of f n a found from Z A r 
determined from a direct fit of C A r (t) to Eq. (|2.6|) . The third column gives f^a found 
from Z A r determined by fitting the ratio C A r {t) /C P r {t) to Eq. (|2.71) and then using 
the value of Z p r found from a fit of C p r (t) to Eq. (|2.4j). The two sets of data in 
all cases are statistically consistent and in all cases, except for the lattice 8 3 x 32 
for k below 0.1500, f n determined from ratio fits has a smaller statistical error f n 
determined from direct fits. The ratio method tends to give less statistical noise, in 
effect, because it uses C p r (t), which is relatively less noisey, to determine m p and 
then extracts only Z A r from the more noisey propagator C A r {t). The direct method 
extracts both Z A r and m p from C A r (t) yielding an m p with greater noise, which 
is then multiplied by a possibly large t and exponentiated, feeding additional noise 
back into the value of Z A r . In the remainder of this paper we use only values of f v 
determined from ratio fits. 

Values of the decay constant in lattice units F p a for the lattices shown in Table 
([!]) are listed in the fourth column of Tables © - (0). The finite renormalizations 
for F p shown in these tables all include both the leading term z vo of Eqs. Q2.ll ) and 
the first order mean- field improved perturbative correction z Y1 of Eqs. ([2.12|) . The 
values of Z^, r used to determine F p were all extracted from direct fits of Cy r (t) to Eq. 
T5T). 
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5 VOLUME DEPENDENCE 



Percentage changes in decay constants going from 8 3 x 32 and 16 3 x 32 to 
24 3 x 32, at j3 of 5.70, are given in Table [11. These changes are the same for all 



choices of finite renormalization. All of the differences appear to be of marginal 
statistical significance and may therefore best be viewed as upper bounds on the 
volume dependence of our results. A variety of different arguments suggest that, 
for the range of k, f3, and lattice volume we have examined, the errors in valence 
approximation decay constants due to calculation in a finite volume L 3 are bounded 
by an expression of the form Ce~ L ^ R . A simple non-relativistic potential model gives 
this expression with the radius of a hadron's wave function for R. A more elaborate 
field theory argument gives for R the Compton wave length of a pair of pions, which 
is the lightest state that can be exchanged between a pair of identical pseudoscalar or 
vector mesons. At (3 of 5.70, R is thus very likely to be between 3 and 5 lattice units. 
Since the changes in decay constants shown in Table [11] between 16 3 and 24 3 are all 
less than 5% for k > 0.1650, it follows that the differences between these values in 
24 3 and in infinite volume should be less than 1.3%. For extrapolations to physical 
quark masses, we use only decay constants with k > 0.1650. 

6 QUARK MASS EXTRAPOLATION 

At the largest k on each lattice, the ratio m^jm p is significantly larger than 
its experimentally observed value of 0.179. Thus to produce decay constants for 
hadrons containing only light quarks, our data has to be extrapolated to larger k 
or, equivalently, to smaller quark mass. We did not calculate directly at larger k 
both because the algorithms we used to find quark propagators became too slow and 
because the statistical errors we found in trial calculations became too large. 

Define the quark mass in lattice units m q a to be 

m « a = k ~ h (61) 

where k c is the critical hopping constant at which m n becomes zero. We found f n a 
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and F p a, for all three possible choices of finite renormalization in Eqs. Q2.10| ) - ( |2.12| ) 
to be nearly linear functions of m q a over the entire range of k considered on each 
lattice. 

Figure [13] shows f n and F p , with finite renormalizations including the first 
order mean- field improved perturbative correction, as functions of m q . Data is shown 
from the three lattices of Table [l] which we use to evaluate continuum limits, 16 3 x 32, 
24 3 x 36 and 30 x 32 2 x 40. For convenience, decay constants at each (3 are shown 
in units of the central value of the rho mass m p (m n ), at the same (3, extrapolated 
to the the "normal" quark mass m n . Quark masses m q for each (3 are shown in 
units of the central value of the strange quark mass m s at the same (3. The data 



Figure [13] has been scaled by the central values of m p {m n ) and m s taken as arbitrary 
external parameters, and the error bars shown do not include the effect of statistical 
fluctuations in m p (m n ) or m s . Table ( |T0| ) gives values of m s found in Ref. |1| by 
requiring m 7r [(m n + m s )/2]/m p (m n ) to be equal to the physical value of mx/rap- 
The lines in Figure [H| are fits of decay constants measured in lattice units, / T a and 
F p a, to linear functions of the quark mass in lattice units, m q a, at the three smallest 
quark masses in the data set at each (3. These fits were obtained by minimizing x 2 
obtained form the full correlation matrix amoung the data points. The correlation 
matrix was calculated by taking averages of data values and products of data values 
over bootstrap ensembles generated as described in Section (|3]). The \ 2 f° r these 
fits, and corresponding fits on the lattice 24 3 x 32 at (3 of 5.70, are given in Table 
|T2| . The fits in Figure [13] appear to be quite good and provide, we believe, a reliable 
method for extrapolating decay constants down to light quark masses. With naive 



finite renormalization, Eq. Q2.10 ), or zero order mean-field finite renormalization, 



Eq. ( j2.ll ), f n a and F p a fit straight lines in m q a about as well as the first order 



perturbatively renormalized data of Figure [L3. 

The linear fits of Figure [13] permit the determination of fx and F^ in addition 
to /„- and F p . For a pion composed of a quark and antiquark with mass m q ^ m q , 



Figure [13] suggests 

f n = a q m q + agm q + (3. (6.2) 
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Charge conjugation invariance then gives a q = a q . It follows that the kaon, which 
is a pion with, say, m q = m s and m q = m n , will have the same decay constant as a 
pion composed of a single type of quark and antiquark with m q = m-q = (m s + m n )/2. 
On the other hand, the linear fits of Figure [13] permit F p to be extrapolated to the 
point m q = m-q = m s which, in the valence approximation, gives F^. Tables |3| and 



14] give the value of f n a, /a- a, F p a and F^a obtained from the fits in Figure [13]. The 
statistical uncertainties in these quantities were obtained by a further application of 
the bootstrap method of Section (|3|). Bootstrap ensembles of the underlying gauge 
configurations were generated, and on each bootstrap ensemble the extrapolated decay 
constants were recalculated. The uncertainty in each decay constant was obtained 
from the resulting distribution. The correlation matrices used to fit bootstrap data 
to linear functions of m q a were taken to be the same as the correlation matrices for 
the full ensemble. To recalulate correlation matrices separately on each bootstrap 
ensemble by a further bootstrap would have been too time consuming. 

7 CONTINUUM LIMIT 

The ratios f n /m p , fx/rap, F p /m p and F^jm p for physical quark masses we 
then extrapolated to zero lattice spacing. For Wilson fermions the leading asymp- 
totic lattice spacing dependence in these decay ratios is expected to be linear in a. 
On the other hand, as shown in Ref. |IJ, m p a follows the two-loop Callan-Symanzik 
scaling prediction in a-ms quite well for the range of (3 considered here. Thus assuming 
the asymptotic form of the lattice spacing dependence of decay ratios appears to be 
reasonable. Figure ffl shows decay constants with first order perturbative renormal- 



ization along with fits to linear functions of m p a. The quantity m p a may be viewed 
as the lattice spacing a measured in units of the physical rho Compton wavelength, 
1/rrip. The vertical bars at m p a of 0, offset slightly for visibility, are the extrapolated 
predictions' uncertainties. The horizontal lines at m p a of lying within or slightly 
above the range of each prediction are the corresponding experimental values. The 
data points in Figure |b| are from the lattices 16 3 x 32, 24 3 x 36 and 30 x 32 2 x 40. 
The values of /3 for these lattices were chosen so that the physical volume in each 
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case is nearly the same. For lattice period L, the quantity m p L is respectively, 9.08 
± 0.13, 9.24 ± 0.19 and, averaged over three directions, 8.67 ± 0.12 |L]. 

The fits shown in Figure [14] were found by minimizing x 2 obtained from the 
full correlation matrix among the fitted data. Since both the x and y coordinates of 
each of the three fitted points on each line have statistical uncertainties, we evaluated 
X 2 among all six pieces of data and chose as fitting parameters the slope and inter- 
cept of the line along with the x coordinate of each point. The required correlation 
matrices were found by the bootstrap method as were the statistical uncertainties of 
the extrapolated predictions. The correlation matrices used in fits for each bootstrap 
ensemble were again taken from the full ensemble and not recalculated on each boot- 
strap ensemble independently. The x 2 P er degree of freedom for the fits in Figure [14], 
are given in Table [15|. 

For the lattice spacing dependence of decay ratios found using zeroth order 
mean-field finite renormalization and using naive finite renormalization we also made 
fits to linear functions of m p a. The slopes with respect to m pa of the ratios f n /m p , 
Jk/^p, F P lm p and F^jm p along with x 2 f° r each fit are also given in Table |H| It is 
clear from these results, as menitioned in Section [|, that naive renormalization leads 
to decay ratios with significantly stronger lattice spacing dependence than found for 
either zeroth or first order improved perturbative renormalization. There also appears 
to be some tendency for first order perturbative renormalization to lead to weaker lat- 
tice spacing dependence than zeroth order. It follows that the extrapolations we have 
done to zero lattice spacing are likely to be most reliable for first order perturbative 
renormalzation and least reliable for naive renormalization. 



8 INFINITE VOLUME LIMIT 



The continuum ratios we found in finite volume were then corrected to infinite 
volume by an adaptation of the method used in Ref. IJ to correct finite volume 
continuum mass ratios to infinite volume. From f^/rrip, with first order perturbative 
renormalization, as a function of lattice spacing a and lattice period L, both measured 
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in physical units, define the finite volume correction term A (a, L) to be 

A(a, L) = — (a, oo) — —(a, L). (8.1) 
m p m p 

The quantity which we would like to determine is A(0, 9/m p ). Now the ratio f^/rrip, 
for L of 9/m.p, undergoes a relative change of a bit less than 20% as a goes from its 
value a 5 7 at (3 of 5.70 to 0. Thus we would expect an error of about 20% of A(0, 9/m p ) 
for the approximation 

A(0, — ) « A(a 5 . 7 , — ). (8.2) 
m p m p 

Moreover, from our earlier discussion of the exponential approach of decay ratios to 
their infinite volume limits, it follows that with an additional error of about 20% of 
A(0,9/m p ) we have 

A(a 5 . 7 , — ) « ^(a 5 . 7 , — ) - A( fl5>7> A), (8 .3) 
m p m p m p m p m p 

where L of 13.5/m p corresponds to the lattice 24 3 x 32 at f3 of 5.70. Finally, a direct 
evaluation of the right side of Eq. (|3.3|) shows it is quite small, with statistical errors 
of about 6% of fn/m p - Combining Eqs. (|8.1|) - ( |8.3|) , we obtain 

A( , oo) * A( , A) + A (Q5 7j 1M) _ A (fl5 7j A). (8 .4) 
m p m p m p m p m p m p m p 

The error in this approximation should be less than about 40% of 6% of f w /m p , 
which is 2.4% of / 7r /m p . Table |l6| lists estimates of the systematic uncertainties in 
equations corresponding to Eq. ( |8.4p for other decay constants and other choices 
of finite renormalization. First order perturbative renormalization consistently gives 
the smallest systematic error in volume correction largely because the lattice spacing 
dependence of these decay ratios is smallest. 

The ratios f n /m p , fx/rrip, F p /m p and F<p/m p , for all three different choices of 
finite renormalization, extrapolated to zero lattice spacing with m p L fixed at 9, and 
then corrected to infinite volume are shown in Table O. For first order perturbative 
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renormalization we also give, in Table [18], finite and infinite volume values of the ratios 
fit I Ik, F p l ' f K and F^j f K . The errors shown for infinite volume ratios are statistical 
only and do not include the estimates we have just given for the systematic error 
in our procedure for making infinite volume corrections. For first order perturbative 



finite renormalization, a comparision of Tables [16] and |17| shows, however, that the 
systematic errors arising from our method of obtaining infinite volume results are 
much smaller than the statistical errors. For the range of /3 used in our extrapolation 
to zero lattice spacing, first order mean-field theory improved perturbation expan- 
sions have been shown |§ to work quite well for a wide variety of different quantities. 
In addition, as we mentioned earlier, the extrapolation to zero lattice spacing should 
be most reliable for this renormalization scheme. Thus we believe that the num- 
bers in Table |17] obtained using first order mean-field perturbative renormalization 
are significantly more reliable than those found using the other two renormalization 
methods. Zeroth order perturbative renormalization has been included, however, to 
provided some measure of the degree to which our results may be sensitive to the 
choice of renormalization. Half of the difference between first order and zeroth order 
mean-field perturbative renormalization appears to us to be a conservative estimate 
of the systematic uncertainty in the first order results arising from the missing sec- 
ond and higher order perturbative renormalization contributions. In all cases this 
uncertainty is significantly less than the statistical errors. Predictions obtained with 



naive renormalization have been included in Table [17] largely as a curiosity. It is 
interesting to notice, however, that the difference between the final, infinite volume 
results found with naive renormalization and those found with first order perturbative 
renormalization is still less than 1.5 times the naive renormalization statistical errors. 



The predicted infinite volume ratios in Table T7 are all statistically consistent 
with the corresponding finite volume ratios. The main consequence of the correction 
to infinite volume is an increase in the size of the statistical uncertainty in each 
prediction. 

The experimental numbers shown in Table [17] for /„- and fx are from charged 
particle decays and for f p from neutral decays. In all cases the uncertainties in 
the experimental values are 0.001 or less. As mentioned in the introduction, an 
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experimental value for the neutral pion decay H gives f n /m p of 0.110 ± 0.005, which 
is quite close to our prediction. The systematic uncertainties in this experimental 
number, however, are larger than those for the charged pion decay. As a result the 
significance of the improved agreement of our prediction with the observed neutral 
pion decay constant is unclear to us. 

We would like to thank Paul Mackenzie for discussions, and Mike Cassera, 
Molly Elliott, Dave George, Chi Chai Huang and Ed Nowicki for their work on GF11. 
We are particularly grateful to Chris Sachrajda for calling our attention to an error 
in an earlier version of this paper. 
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lattice 


P 


k 




skip 


count 


8 3 x 32 


5.70 


0.1400 - 


0.1650 


1000 


2439 


16 3 x 32 


5.70 


0.1600 - 


0.1675 


2000 


219 


24 3 x 32 


5.70 


0.1600 - 


0.1675 


4000 


58 


24 3 x 36 


5.93 


0.1543 - 


0.1581 


4000 


210 


30 x 32 2 x 40 


6.17 


0.1500 - 


0.1532 


6000 


219 



Table 1: Configurations analyzed. 



lattice 


P 


L/ 02 




O 02/°02 


°02 


16 3 x 32 


5.70 


0.09 


0.22 


0.23 


0.11 


24 3 x 36 


5.93 


0.08 


0.07 


0.47 


0.68 


30 x 32 2 x 40 


6.17 


0.16 


0.10 


0.01 


0.25 



Table 2: Values of x 2 per degree of freedom for fits to propagators. 



lattice (3 k c 



8 3 


x 32 


5.70 


0.169012 ±0.000102 


16 3 


x 32 


5.70 


0.169405 ±0.000052 


24 3 


x 32 


5.70 


0.169304 ±0.000035 


24 3 


x 36 


5.93 


0.158948 ±0.000026 


30 x 32 2 


x 40 


6.17 


0.153763 ±0.000018 



Table 3: Critical hopping constant values. 
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p 


Z A1 


z vi 


5.70 


0.93269 


0.82194 


5.93 


0.94208 


0.84679 


6.17 


0.94828 


0.86320 



Table 4: First order perturbative correction factors to finite renormalizations. 



k f n a direct f n a ratio F p a 

0.1400 0.2346 ±0.0011 0.2347 ± 0.0011 0.2861 ± 0.0015 

0.1450 0.2080 ±0.0014 0.2083 ± 0.0010 0.2684 ± 0.0014 

0.1500 0.1811 ±0.0016 0.1820 ±0.0010 0.2505 ± 0.0017 

0.1550 0.1532 ±0.0018 0.1538 ± 0.0012 0.2316 ± 0.0026 

0.1600 0.1238 ±0.0022 0.1243 ± 0.0012 0.2074 ± 0.0022 

0.1650 0.0895 ± 0.0074 0.0925 ± 0.0021 0.1843 ± 0.0043 

Table 5: Decay constants in lattice units measured on a lattice 8 3 x 32 at (3 = 5. 
Finite renormalizations include the first order perturbative correction. 
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k f w a direct f n a ratio F p a 

0.1600 0.1212 ±0.0040 0.1236 ± 0.0016 0.1975 ± 0.0047 

0.1650 0.0915 ±0.0046 0.0947 ± 0.0019 0.1840 ± 0.0069 

0.16625 0.0799 ±0.0061 0.0875 ± 0.0022 0.1794 ± 0.0061 

0.1675 0.0720 ± 0.0098 0.0818 ± 0.0025 0.1794 ± 0.0034 

Table 6: Decay constants in lattice units measured on a lattice 16 3 x 32 at (3 = 5.70. 
Finite renormalizations include the first order perturbative correction. 



k f n a direct f n a ratio F p a 

0.1600 0.1251 ±0.0026 0.1249 ± 0.0010 0.2116 ± 0.0035 

0.1650 0.0951 ±0.0031 0.0940 ± 0.0010 0.1930 ± 0.0055 

0.1663 0.0792 ± 0.0083 0.0860 ± 0.0015 0.1865 ± 0.0052 

0.1675 0.0769 ± 0.0059 0.0787 ± 0.0018 0.1752 ± 0.0043 

Table 7: Decay constants in lattice units measured on a lattice 24 3 x 32 at (3 = 5.70. 
Finite renormalizations include the first order perturbative correction. 



k f w a direct f T a ratio F p a 

0.1543 0.0763 ± 0.0025 0.0729 ± 0.0009 0.1187 ± 0.0030 

0.1560 0.0668 ±0.0037 0.0625 ± 0.0008 0.1096 ± 0.0043 

0.1573 0.0562 ±0.0021 0.0544 ± 0.0009 0.1077 ± 0.0033 

0.1581 0.0517 ±0.0029 0.0490 ± 0.0012 0.1027 ± 0.0037 

Table 8: Decay constants in lattice units measured on a lattice 24 3 x 36 at (3 = 5.93. 
Finite renormalizations include the first order perturbative correction. 
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k f w a direct f n a ratio F p a 

0.1500 0.0557 ±0.0009 0.0566 ± 0.0005 0.0852 ± 0.0008 

0.1519 0.0444 ±0.0012 0.0449 ± 0.0004 0.0760 ± 0.0010 

0.1526 0.0403 ±0.0017 0.0400 ± 0.0005 0.0724 ± 0.0013 

0.1532 0.0361 ±0.0023 0.0356 ± 0.0007 0.0696 ± 0.0028 

Table 9: Decay constants in lattice units measured on a lattice 30 x 32 2 x 40 at 
j3 = 6.17. Finite renormalizations include the first order perturbative correction. 



lattice j3 m p a m n a m s a 

16 3 x 32 5.70 0.5676 ±0.0079 0.00390 ± 0.00012 0.09662 ± 0.00291 

24 3 x 32 5.70 0.5409 ± 0.0089 0.00348 ± 0.00012 0.08620 ± 0.00303 

24 3 x 36 5.93 0.3851 ± 0.0079 0.00223 ± 0.00009 0.05536 ± 0.00226 

30 x 32 2 x 40 6.17 0.2768 ± 0.0039 0.00141 ± 0.00004 0.03503 ± 0.00096 

Table 10: The normal and strange quark masses, m n and m s , and the rho mass 
extrapolated to quark mass m n . 
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decay k /(24 3 ) - /(8 3 ) /(24 3 ) - /(16 3 ) 

U 0.1600 ~ 0.4 ±1.4% ~ 1.0 ±1.5% 

0.1650 1.7 ±2.6% -0.7 ±2.3% 

0.1663 -1.7 ±3.3% 

0.1675 -4.1 ±4.1% 

' _ ¥ p 0.1600 2.0 ±1.7% 6.7 ±2.9% 

0.1650 4.5 ±3.9% 4.7 ±5.1% 

0.1663 3.8 ±4.6% 

0.1675 -2.4 ±2.8% 

Table 11: Changes in decay constants from 8 3 x 32 to 24 3 x 32 and from 16 3 x 32 to 
24 3 x 32 at (3 = 5.70 



lattice (3 \ 2 f° r fw X 2 f° r F p 

16 3 x 32 5.70 0.256 0.149 

24 3 x 32 5.70 0.003 0.180 

24 3 x 36 5.93 0.394 0.358 

30 x 32 2 x 40 6.17 0.130 0.014 

Table 12: Values of \ 2 f° r fits of decay constants to linear functions of m q at the three 
smallest values of m q . 
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lattice (3 



16 3 x 32 5.70 0.0724 ±0.0030 0.0861 ± 0.0017 

24 3 x 32 5.70 0.0691 ± 0.0025 0.0831 ± 0.0016 

24 3 x 36 5.93 0.0449 ± 0.0012 0.0532 ± 0.0008 

30 x 32 2 x 40 6.17 0.0323 ± 0.0008 0.0378 ± 0.0005 

Table 13: Decay constants in lattice units extrapolated to physical quark mass. 



lattice (3 F p a F^a 

16 3 x 32 5.70 0.1747 ±0.0076 0.1870 ± 0.0107 

24 3 x 32 5.70 0.1642 ± 0.0070 0.1975 ± 0.0063 

24 3 x 36 5.93 0.1014 ± 0.0047 0.1098 ± 0.0037 

30 x 32 2 x 40 6.17 0.0668 ± 0.0026 0.0748 ± 0.0012 

Table 14: Decay constants in lattice units extrapolated to physical quark mass. 
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decay 


renorm. 


slope 


x 2 




first order perturb. 


0.035 ±0.022 


0.62 




zeroth order perturb. 


0.044 ±0.023 


0.61 




naive 


0.112 ±0.030 


0.74 


fK/m p 


first order perturb. 


0.052 ±0.016 


1.04 




zeroth order perturb. 


0.063 ±0.017 


1.02 




naive 


0.112 ±0.022 


0.99 


Fp/m p 


first order perturb. 


0.228 ± 0.047 


0.04 




zeroth order perturb. 


0.325 ±0.058 


0.06 




naive 


0.546 ± 0.075 


0.19 


F^/nip 


first order perturb. 


0.192 ±0.061 


0.20 




zeroth order perturb. 


0.286 ±0.073 


0.22 




naive 


0.362 ±0.091 


0.20 



Table 15: Slope and \ 2 of fits of decay ratios' lattice spacing dependence to linear 
functions of m p a. 
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decay renorm. 



error 



U/m p 


first order perturb. 


2% 




zeroth order perturb. 


3% 




naive 


4% 


fK/m p 


first order perturb. 


2% 




zeroth order perturb. 


2% 




naive 


2% 


Fp/m p 


first order perturb. 


6% 




zeroth order perturb. 


7% 




naive 


10% 




first order perturb. 


8% 




zeroth order perturb. 


10% 




naive 


12% 



Table 16: Estimated systematic errors in infinite volume quantities arising from the 
correction from finite volume to infinite volume. 
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decav 


renorm. 


finite volume 


infinite volume 


obs. 


f TYI 
Jit/ " L p 


first order nertnrb 


106 ± 009 


106 ± 014 


0.121 




7,prnt,b nrHpr nprtnrh 


110 ± 009 


110 ± 015 






naive 


0.119 ± 0.011 


0.120 ± 0.019 




fK/m p 


first order perturb. 


0.121 ±0.006 


0.123 ±0.009 


0.148 




zeroth order perturb. 


0.125 ±0.007 


0.127 ±0.010 






naive 


0.141 ±0.008 


0.144 ±0.013 




Fp/m p 


first order perturb. 


0.177 ±0.021 


0.173 ±0.029 


0.199 




zeroth order perturb. 


0.189 ±0.024 


0.184 ±0.036 






naive 


0.191 ±0.030 


0.186 ±0.045 




F^/rrip 


first order perturb. 


0.217 ±0.019 


0.253 ±0.035 


0.219 




zeroth order perturb. 


0.234 ±0.023 


0.277 ±0.043 






naive 


0.269 ±0.029 


0.328 ± 0.052 





Table 17: Calculated values of meson decay constants extrapolated to zero lattice 
spacing in finite volume, then corrected to infinite volume, compared with observed 
values. 



decay 


finite volume 


infinite volume 


obs. 


/tt/ Jk 


0.875 ± 0.034 


0.864 ±0.065 


0.818 


F p /fK 


1.467 ±0.170 


1.412 ±0.229 


1.345 


f*ISk 


1.795 ±0.131 


2.058 ±0.283 


1.480 



Table 18: Calculated values of meson decay constants, using first order perturbative 
renormalization, extrapolated to zero lattice spacing in finite volume, then corrected 
to infinite volume, compared with observed values. 
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Figure 1: Effective masses and fitted mass for the pseudscalar propagator C^ 2 (t) on 
the lattice 16 3 x 32 at at (3 = 5.70 and k = 0.1675 
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Figure 2: Effective masses and fitted mass for the axial vector propagator Q^(i) on 
the lattice 16 3 x 32 at at (5 = 5.70 and k = 0.1675 
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Figure 3: The propagator ratio C^t) j 'C^(i) and a fit to its large t plateau on the 
lattice 16 3 x 32 at at (5 = 5.70 and k = 0.1675 
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Figure 4: Effective masses and fitted mass for the vector propagator C^ 2 (t) on the 
lattice 16 3 x 32 at at (5 = 5.70 and k = 0.1675 
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Figure 5: Effective masses and fitted mass for the pseudscalar propagator C^ 2 (t) on 
the lattice 24 3 x 36 at at (3 = 5.93 and k = 0.1581 
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Figure 6: Effective masses and fitted mass for the axial vector propagator Q^(i) on 
the lattice 24 3 x 36 at at (3 = 5.93 and k = 0.1581 
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Figure 7: The propagator ratio C^(t) / C^(i) and a fit to its large t plateau on the 
lattice 24 3 x 36 at at f3 = 5.93 and k = 0.1581 
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Figure 8: Effective masses and fitted mass for the vector propagator C^ 2 (t) on the 
lattice 24 3 x 36 at at (3 = 5.93 and k = 0.1581 
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Figure 9: Effective masses and fitted mass for the pseudscalar propagator C^(i) on 
the lattice 30 x 32 2 x 40 at at (5 = 6.17 and k = 0.1532 
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Figure 10: Effective masses and fitted mass for the axial vector propagator C^ 2 (t) 
on the lattice 30 x 32 2 x 40 at at (5 = 6.17 and k = 0.1532 
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Figure 11: The propagator ratio C^it) /C^if) an d a fit to its large t plateau on the 
lattice 30 x 32 2 x 40 at at f3 = 6.17 and k = 0.1532 
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Figure 12: Effective masses and fitted mass for the vector propagator C^ 2 (t) on the 
lattice 30 x 32 2 x 40 at at (3 = 6.17 and k = 0.1532 
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Figure 13: The decay constants F p and with perturbative renormalization, in 
units of the central value of the physical rho mass m p (m n ), as functions of the quark 
mass m q , in units of the central value of the strange quark mass m s . 
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Figure 14: Perturbatively renormalized decay constants as functions of the lattice 
spacing a, in units of l/m p . The error bars near zero lattice spacing are uncertainties 
in the extrapolated ratios, and the horizontal lines represent experimentally observed 
values. 



40 



